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(57) Abstract 

A radio signal travels over a radio channel 
in the form of muUipath propagation, whereby a 
signal received by the base station and transmit- 
ted by the terminal equipment is composed of 
multipath components propagated over different 
distances and thus reaching the receiver at dif- 
ferent times. Thus, the information required for 
computing the bearing of the received signal can 
be extracted from the arrival time differences of 
the multipath propagated signal- An antenna ar- 
ray formed by a plurality of antenna elements 
receives first the multipath propagated signal. 
Next, means (41. 4N) for determination of 
channel impulse response form an estimate of 
the radio channel impulse response from the re- 
ceived signal of each antenna array element, af- 
ter which correlation means (46) correlate each 
channel impulse response estimate with a refer- 
ence estimate selected from the group of said 
channel impulse response estimates. Further, 
said means choose as the time difference of the 
received signal that delay value giving the maxi- 
mum correlation. Finally, computational means 
(46) estimate the bearing of the radio signal 
transmitted by the terminal equipment from the 
arrival time differences thus determined. Using 
a two-dimensional antenna array, both the bear- 
ing sector and elevation of the received signal 
can be estimated. 
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DETERMINATION OF TERMINAL LOCATION IN A RADIO SYSTEM 

Scope of the invention . 

The present invention relates to a method of locating the position 
5 of terminal equipment in a system in which one part of a subscriber connec- 
tion is formed by a radio connection between the terminal equipment and a 
base station. Such a system can be, e.g., a cellular mobile radio system of 
several base stations in which the terminal equipment during a connection 
may move from one cell to another. 

10 

Background of the invention 

Conventionally, information on the accurate position of terminal 
equipment during the connection between the terminal equipment and a 
base station has not generally been considered crucially important. However, 
15 in planning of radio systems such as cellular telephone systems, the dis- 
tance of the terminal equipment from the base station must be taken into 
account, with particular attention being paid on the variation of the distance 
during connection. This is vital especially in time division digital cellular sys- 
tems in which transmission is carried out in bursts during a given transmit 
20 time slot. In order to prevent bursts transmitted at the same carrier frequency 
in consecutive time slots from over-lapping in the received signal at the base 
station, the bursts transmitted by each individual mobile terminal are ex- 
tended by a guard time of a certain duration at the end of each burst. Such a 
guard time is necessary, because the mobile terminals transmitting at the 
25 same radio carrier, but in consecutive time slots, are located at random dis- 
tances from the base station, whereby also the propagation time of radio 
waves from the base station to the mobile terminal will vary from time slot to 
time slot. The guard time is determined so that the base station, within the 
timing margin of the transmit/receive frame, measures the arrival time differ- 
30 ence between the burst transmitted by the base station and the burst re- 
ceived from the mobile terminal and subsequently computes a timing ad- 
vance (TA), after which the base station transmits the value thereof to the 
mobile terminal. Thus, the base station dynamically controls the transmit time 
of each mobile terminal on the basis of bursts received therefrom. In a GSM 
35 system, the timing advance may be assigned values in the range 0-233 ms, 
which in relation to the time slot overall duration means that a mobile termi- 
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nal operating with tine maximum timing advance cannot be fartlier away than 
35 km from the base station. 

As it is generally known, the location of an emitting signal source 
always requires more than one receiver, spatially separated from each other. 
5 Methods of signal direction estimation are based on the estimation of time 
differentials needed by the same signal to reach the different receivers. Es- 
timation of received signal arrival time differences is performed by correlating 
the output signals from different receivers. The intersignal delay which gives 
the maximum value of the correlation function represents the estimate of 
10 maximum likelihood for the arrive time difference, provided that said value of 
the correlation function exceeds a preset threshold value. 

In a cellular radio system, the above-described procedure is car- 
ried out so that the transmit signal of a mobile terminal is received at a plu- 
rality of base stations, or alternatively, at a base station equipped with an 
15 array of antenna elements with individual receivers. Then, the received sig- 
nals obtained from the different base stations or an antenna array of a single 
base station can be correlated with each other in order to determine the 
arrival time differences, and therefrom, the bearing angle of the received 
signal. 

20 These conventional approaches are hampered by the computa- 

tional complexity involved in relation to the accuracy offered by them. 

However, in some cases it is a desirable property of the network 
to provide maximally exact information on the location of the mobile terminal. 
One such case is encountered when it is required to restrict the movement of 
25 the mobile terminal to occur within the coverage area of a given cell only or 
when a more detailed scheme of call charging rates is desired. 

One approach to the location of a mobile terminal is presented in 
Fl Pat. Appl. No. 963,382, filed in Aug. 30, 1996, not yet public at the filing 
date of the present application. This method uses an estimated value of the 
30 timing advance, which is a readily available value computed by the commu- 
nications system, proportional to the distance between the terminal equip- 
ment and the base station. \A/hile the estimate of the timing advance com- 
puted for the distance between the serving base station and the terminal 
equipment gives only a rough estimate of the terminal equipment's distance 
35 from the base station, but the estimate of the terminal equipment timing 
advance computed for a greater number of base stations already gives the 
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terminal equipment location with a relatively good accuracy. The latter tech- 
nique is supported by the OTD (Observed Time Difference) property of 
Phase 2 of the GSM standard that directly gives the propagation time differ- 
ences of the signals from the mobile terminal to the different base stations. 
5 Cited patent application further suggests that, in CDMA (Code Division Multi- 
ple Access) systems, the mobile terminal location can be determined by 
correlating the pilot signals transmitted by the base stations. 

The location of terminal equipment is especially important also in 
authority networks. These networks are private cellular networks which are 

10 particularly designed to operate under different distress situations, too. 
Herein, it is important that the network can identify the location of the termi- 
nal equipment without the need for the mobile terminal operator to verbally 
report his location. In open networks, the mobile terminal location information 
would also be a desirable property in emergency call situations, because the 

15 emergency call originator does not necessarily always know his exact loca- 
tion. 

It is an objective of the present invention to provide a mobile ter- 
minal location method which is characterized by simple computation routines 
and maximum utilization degree of the properties of existing communication 
20 networks. 

The goals of the invention are achieved by virtue of the specifica- 
tions expressed in the appended independent claims. 

Summary of the invention 

25 The invention is based on the idea that, since a radio signal trav- 

els over a radio channel in the form of multipath propagation, whereby the 
terminal equipment transmit signal received by the base station is composed 
of time-dispersed multipath signal components propagated over different 
distances and thus reaching the receiver at different times, all the information 

30 required for computing the bearing of the received signal can be extracted 
from the differences of the dispersed arrival times of the multipath signal. 
Accordingly, the multipath propagated signal is first received by means of an 
antenna array of the base station, whereafter the channel estimates of the 
signals received by the individual antenna array elements are formed and 

35 then correlated to determine the relative time differences between them. 
Finally, the estimated bearing of the received signal can be computed from 
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these differences of the arrival time estimates of the individual antenna sig- 
nals. 

The method according to the invention is suited for use in con- 
junction with a conventional one-dimensional antenna array, so called a 
5 linear array, whereby only the bearing angle of the received signal can be 
estimated. According to a preferred embodiment of the invention, a two- 
dimensional or three-dimensional antenna array is used. Then, the bearing 
angle information can be complemented with the estimation of the signal 
elevation angle. Hence, using geometrical rules, the location of the terminal 
10 equipment can be computed from the known values of the antenna array 
height from the ground and the estimated bearing and elevation angles of 
the signal received by the antenna array. 

The accuracy of position determination may further be improved if 
the concept can be extended to a greater number of base stations, each 
15 equipped with a same type of antenna array. Then, the estimate of received 
signal bearing at each base station is determined according to the invention, 
and the actual location of the mobile terminal is subsequently computed from 
these bearing angle estimates. 

A significant benefit of the disclosed method over conventional 
20 techniques is its computational simplicity. As only the channel estimates are 
correlated, the number of required computational operations is reduced radi- 
cally. The resolution achieved herein is at least as good as with conventional 
methods, because the arrival time differences between the individual an- 
tenna signals can be determined with the resolution of one sample interval, 
25 and the lesser complexity of arrival time difference determination can be 
utilized in the postprocessing of the information. Particularly with the use of 
two or three-dimensional antenna arrays, postprocessing can render a sig- 
nificant improvement in the achievable accuracy, because the bearing of the 
received signal changes relatively slowly. 

30 

Brief description of drawings 

In the following, the invention will be examined in greater detail 
with reference to the accompanying drawings, in which 

35 Figure 1 depicts a cellular network; 

Figure 2 shows an antenna array at a base station; 
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Figure 3 shows a graph of a delay distribution of the received signal; 

Figure 4 shows the block diagram of a receiver; 

Figure 5 shows some channel response estimates; 

Figure 6 depicts correlation of channel response estimates; 

Figure 7 show$ a flow chart of the method; and 

Figure 8 illustrates determination of mobile terminal distance. 



Detailed description of the invention 

Referring to Fig. 1, therein is illustrated a typical cellular system, in 
which the terminal equipment, here a mobile station MS, is connected via a 
radio channel to a base station BTS1, from where the connection is contin- 
ued further over wire. The terminal equipment may move during the connec- 
tion from one cell to Another, requiring the radio channel to be correspond- 
ingly handed over from one base station to another, e.g., in a sequence 
BTS1^BTS2-^BTS3. 

A method conventionally used in a TDMA system for locating ter- 
minal equipment is based on the use of estimated timing advance values 
TA1, TA2 and TA3. As the value of each timing advance is proportional to 
the terminal equipment distance from the base station, it is easy to estimate 
the terminal equipment location relatively accurately on the basis of these 
timing advance estimates. This approach has a particular advantage in that 
the communication system anyhow computes the timing advance values. 

The present invention draws upon the idea of using a conven- 
tional antenna array. Such an antenna array is comprised of a plurality of 
antenna elements acting as individual antennas. By phasing the antenna 
elements in this type of antenna array, it is possible to steer the maximum of 
the radiation power in a desired direction or spatial angle, while in other 
directions the radiation power is strongly suppressed, i.e the antenna array is 
used in forming a desired radiation pattern. 

The simplest configuration of an antenna array is the one- 
dimensional linear array having the elements equidistantly spaced in a linear 
fashion; the next more complex configuration is a two-dimensional array 
having the elements located in the same plane, e.g., at the corners of a 
rectangle; and the most complex configuration is a three-dimensional array 
having the elements forming a spatial structure such as a cube or pyramid. 
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Arrangements are provided for individual handling of the element re- 
ceive/transmit signals, respectively, in each antenna configuration. 

Now referring to Fig. 2, an exemplifying embodiment of antenna 
array mounted at the top of an antenna mast 22 has a three-dimensional 
5 configuration comprising a first set of antenna elements a, b, c, d, e, f, k, I 
in an xy plane at the intersection points of an equidistantly spaced grid, and 
a corresponding second set of antenna elements m, n, etc., also in an xy 
plane, spaced vertically at a distance z from the first plane. The signals re- 
ceived and transmitted by the elements are processed in an antenna signal 

10 processing block 21. 

In the description below, it is assumed that the cellular system 
concerned is a CDMA system, whereby the base stations of Fig. 2 perform 
the transmission by spreading the information over the entire transmission 
bandwidth by means of a spreading code individually allocated for each 

15 connection. Thus, the terminal equipment such as a mobile station MS (Fig. 
1), for instance, has established a connection with a base station BTS1 by 
spreading its transmission power over the entire radio channel bandwidth 
using its individual spreading code. Additionally, it is assumed that the base 
station is equipped with a two- or three-dimensional antenna array, e.g., of 

20 the type shown in Fig. 2. 

The CDMA signal transmitted by the terminal equipment, whose 
location is to be determined, reaches each element of the antenna array in 
the base station as a plane wave. The signal received by each antenna array 
element has undergone multipath propagation, thus beirig composed of a 

25 number of time dispersed signal components reaching each element at dif- 
ferent times. From each antenna element, the signals are further taken to a 
conventional Rake receiver. A separate receiver is provided for each an- 
tenna element. 

in Fig. 4 is shown a receiver comprising a bank of N correlator 
30 blocks, denoted in the diagram as blocks 41, 42, 4N. In each block, the 
received signal taken to the block from an antenna element is correlated with 
the same code used for spreading the signal in the sending terminal equip- 
ment, but by varying the time delay, whereby the multipath propagation 
components can be obtained from the received signal. According to the 
35 number N of the correlator blocks, N pes. of multipath propagation compo- 
nents can be resolved that in combination form the impulse response of the 
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radio channel between a given antenna array element and the terminal 
equipment. 

The channel impulse response may have a shape similar to the 
one shown in Fig. 3 in the case that 4 pes. correlator blocks (N=4) are used. 
5 Here, tap 1 is produced by correlator block 1 , tap 1 by correlator block 2, etc. 
The heights of the taps correspond to the spectral power of each signal 
component propagated over different paths, and the distances of the taps 
are proportional to the arrival time difference of the impulse signal compo- 
nent relative to a zero time point representing a reference time. Accordingly, 
10 the arrival delay to tap 1 is DT1, the arrival delay to tap 2 is PT2, the arrival 
delay to tap 3 is DT3, and the arrival delay to tap 4 is DT4. Corresponding 
multipath channel response estimates are obtained from each receiver con- 
nected to an antenna array element and all the signal estimates are taken to 
a combinatory block, where they are summed coherently. The summing 
15 result is then passed to a detector. 

For the most of its parts, the above-described technique is well 
known in the art. Now, according to the invention, the conventional summing 
step is complemented by correlating the multipath channel response esti- 
mates produced by each of the receivers. The multipath channel response 
20 estimates and their propagation delay values obtained from the correlator 
bank are taken to a correlator/position computing block 46. As the estimation 
of the received signal bearing is based on the estimation of the arrival time 
differences between the instants when the same signal reaches the different 
receivers, the estimation of these time differences can be performed in the 
25 block 56 by correlating the output signals of the different receivers. Hence, 
the input signals of the block 56 also include the channel response estimates 
from the outputs of the other receivers of the antenna array elements. The 
delay value between the channel response estimates producing the highest 
correlation gives an estimate for the arrival time difference, provided that said 
30 correlation function value is larger than a preset threshold. 

Next, the function of the invention is explained by making refer- 
ence to Fig. 7. Correlation is advantageously carried out so that a multipath 
channel response estimate Eref obtained from any one of the antenna array 
elements is selected as the reference, stage 71. For instance, in the diagram 
35 of Fig. 2, the multipath channel response processed from the antenna signal 
of array element I has been selected as the reference. Next in the procedure, 
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a search routine is carried out to test how much the multipath channel re- 
sponse estimate obtained from the antenna signal of the adjacent array 
element k nnust be shifted temporally to achieve the best coincidence 
(matching) with the multipath channel response estimate Ep^p obtained from 
5 the antenna signal of array element I, stage 72. After this correlation is 
maximized, the corresponding value of the temporal shift tj, which is equiva- 
lent to the arrival delay difference, is stored, stage 73, and the equivalent 
correlation operation is performed for the next antenna array element j in 
sequence. Thus, the correlation routine is sequentially repeated for each 
10 multipath channel response estimate produced by the correlator block, after 
which the best match of the multipath channel response estimate for the 
received signal of each antenna array element is obtained. It must be noted 
herein that the receivers connected to the antenna array elements must 
operate in a synchronous mode in order to obtain a correct result. 
15 After each of the antenna array element received signals are 

processed, it is known how much the channel response estimate of each 
antenna array element must be temporally shifted with respect to the tap 
diagram (refer to the tap diagram of Fig. 3) of the reference element, 
vvhereby the delay difference of the received signal in each antenna array 
20 element with respect to the reference element received signal is resolved. 
Subsequently, the bearing of the antenna signal can be computed, stage 75. 

The correlation procedure is further elucidated in Fig, 5, whose 
upper part shows the four-tap diagram of the channel response estimate for 
the received signal of reference element 1, while the lower part shows the 
25 corresponding four-tap channel response estimate for antenna array element 
e. The channel response estimates are plotted on the same time scale, and 
the common time reference point is formed by the reference time point of 
channel response signal taps of the reference element itself. The function of 
correlation is to find such a value of time shift which, when applied to the 
30 channel response estimate of element e, maximizes the value of each de- 
sired correlation function. 

In Fig. 6 is illustrated a possible shape of the correlation function 
when the channel response estimate of antenna array element e is shifted 
along the time axis. It can be seen from the graph that the best correlation is 
35 obtained when the channel response estimate of element e is shifted by a 
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time te- Obviously, this time is also equivalent to the arrival time difference 
between the received signals at the antenna array elements. 

After the arrival time differences between each individual antenna 
element received signal is computed, the bearing of the received signal can 
5 be estimated using the geometry notation of Fig. 8. If the antenna array used 
is a one-dimensional linear array, only the sector of the received signal, that 
is, its bearing angle a can be determined, while the distance of the terminal 
equipment cannot be computed. In most cases, information on the bearing 
sector of the terminal equipment will suffice. On the other hand, if the an- 
10 tenna array is two- or three-dimensional, also the signal elevation angle b of 
the received signal can be computed. Thus, on the basis of the known spa- 
tial angle coordinates (bearing and elevation) of the signal source, normal 
rules of geometry can be applied for computing the distance d of the terminal 
equipment from the base station. Accordingly, the procedure is capable of 
15 locating the terminal equipment. 

As described above, the present method is based on extracting all 
the information required for the determination of the received signal bearing 
from the arrival time differences of the multipath signals, whereby other in- 
formation carried by the received signal, such as the subscriber transmit 
20 signal spreading code or the transmitted symbol data, does not essentially 
affect the estimation of the received signal bearing. 

The most advantageous application of the invention is found in a 
CDMA system, because the multipath channel response estimates are con- 
tinuously and directly available from the CDMA receiver, thus obviating a 
25 separate computation thereof. Another possible application can be found in a 
TDMA system in which the channel impulse response of the received signal 
is computed for the purpose of adjusting the channel equalizer. The multi- 
element antenna array provides-a number of channel responses which are 
correlated with respect to a common time reference point, after which the 
30 multipath propagation time differences can be computed. Finally, the bearing 
of the received signal can be computed from the arrival time differences. 

In both applications, the computational accuracy can be improved 
by subjecting the results to processing such as low-pass filtration, for in- 
stance, since the received signal bearing under realistic conditions varies at 
35 a relatively slow rate. Thus, the filtration process permits elimination of ran- 
dom transients from the estimates of the mobile station distance. 
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Analogously with all other, techniques known in the art, the accu- 
racy of the present nnethod is largely dependent on the quality of the radio 
channel. If a direct line-of-sight can be established between the base station 
and the terminal equipment, the location of the mobile station is easy to 
5 determine, but under conditions of large channel response dispersion, 
whereby the received signal is comprised of multi reflected signal compo- 
nents, the location of the terminal equipment may become an impossible 
task. Such a situation may occur in large macrocells of a network. 

The above description of the method according to the invention 

10 concerns an application to a single base station only. When desired, the 
accuracy of terminal equipment location can be improved by exploiting a 
technique called macrodiversity, which means an extension of the computa- 
tional grid to cover the neighbouring base stations, too. Thus, the mobile 
terminal location can be performed using the bearing estimates computed by 

15 each of the base stations. This arrangement presumes that the bearing in- 
formation is gathered to a single point, where the final computation of the 
terminal equipment location is performed. In this arrangement, the receivers 
of the different base stations need not necessarily be run in mutual synchro- 
nism, but the base stations should have identical antenna arrays and the 

20 receivers connected to the antenna array must be run in mutual synchronism 
at each base station. 
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Claims 

1. A method of locating the position of terminal equipment in a 
system in which one part of the subscriber connection is formed by a radio 
connection between the terminal equipment and a base station and in which 

5 a signal transmitted by the terminal equipment is received in the form of a 
multipath propagated signal by the base station antenna, said antenna com- 
prising an antenna array formed by a plurality of antenna elements, 
characterized in that 

an estimate of the radio channel impulse response is formed from 
10 the received signal of each antenna array element, 

each channel impulse response estimate is correlated with a ref- 
erence estimate selected from the group of said channel impulse response 
estimates and the propagation time difference of each channel impulse re- 
sponse component is set equal to the propagation delay value giving the 
15 maximum correlation, and 

the bearing of the radio signal transmitted by the terminal equip- 
ment is estimated from the time differences thus determined. 

2. A method as defined in claim 1, characterized in that 
the antenna array used is an at least two-dimensional array, whereby the 

20 direction of the radio signal can be resolved for both the bearing sector and 
elevation angle of the received signal. 

3. A method as defined in claim 1 , characterized in that 
the antenna array used is a one-dimensional array, whereby the direction of 
the radio signal can be resolved only for the bearing sector of the received 

25 signal. 

4. A method as defined in claim 1, characterized in that 
said system is a CDMA system and that the channel impulse response esti- 
mates directly available from the CDMA receiver are utilized. 

5. A method as defined in claim 1, characterized in that 
30 the radio signal transmitted by said terminal equipment is received 

at a plurality of base stations, and 

the location of said terminal equipment is resolved on the basis of 
bearing information estimated at each base station. 

6. A method as defined in claim lorS, characterized in 
35 that the resolved bearing information is subjected to low-pass filtration. 
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7. A system in which a part of the subscriber connection is formed 
by a radio connection between the terminal equipment and a base station 
and in which system a signal transmitted by the terminal equipment is re- 
ceived in the form of a multipath propagated signal by the base station an- 

5 tenna, which comprises an antenna array formed by a plurality of antenna 
elements, 

characterized in that said system further includes 
means (41, 4N) for determination of channel impulse response, 
said means forming an estimate of the radio channel impulse response from 
10 the received signal of each antenna array element, 

- correlations means (46), said means correlating each channel 
impulse response estimate with a reference estimate selected from the group 
of said channel impulse response estimates and, for each antenna array 
element, said correlation means choosing as the time difference of the re- 

15 ceived signal that delay value giving the maximum correlation, and 

- computation means (46) for estimating the bearing of the radio 
signal transmitted by the terminal equipment from the arrival time differences 
thus determined. 

8. A system as defined in claim 7, characterized in that 
20 said antenna array is an at least two-dimensional array, whereby the direc- 
tion of the radio signal can be resolved for both the bearing sector and ele- 
vation of the received signal. 

9. A system as defined in claim 7, characterized in that 
the antenna array used is a one-dimensional array, whereby the direction of 

25 the radio signal can be resolved only for the bearing sector of the received 
signal. 

1 0. A method as defined in claim 7, characterized in that 
said system is a CDMA system and that said means (41, 4N) for deter- 
mining the channel impulse response are an inherent part of a conventional 

30 CDMA receiver. 
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